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Abstract 

The probabiUty that an interval I is free of eigenvalues in a matrix ensemble with unitary 
symmetry is given by a Fredholm determinant. When the weight function in the matrix ensemble is 
a classical weight function, and the interval / includes an endpoint of the support, Tracy and Widom 
have given a formalism which gives coupled differential equations for the required probability and 
some auxilary quantities. We summarize and extend earlier work by expressing the probability 
and some of the auxilary quantities in terms of Painleve transcendents. 

Dedicated to Martin Kruskal on the occasion of his 75th Birthday. 



1 Motivations from Physics and Mathematics 

Diverse branches of physics and mathematics provide many examples of quantities which are known to 
have the same statistical properties as the eigenvalues of large random matrices. This is the case for 
quantum chaos — the study of the quantum mechanics of systems which exhibit chaos at the classical 
level. Here random matrices have emerged as a paradigm, in that those systems for which the underlying 
classical mechanics is chaotic are observed to have energy levels with the same statistical properties 
as the eigenvalues of large random matrices, while systems with integrable classical dynamics do not. 
Many examples are known from the study of two-dimensional billiard systems — a single particle in a 
closed domain in two dimensions. One of the most celebrated is the Sinai billiard, which is a square 
billiard with a hard wall circular scatterer in the centre. Numerical studies show that this system, which 
is known to be chaotic at the classical level, has energy levels with the same statistical properties of 
large real symmetric matrices (see e.g. Q). The status of our understanding is reflected in a conjecture 
by Bohigas-Giannoni-Schmit which states: 

"The spectra of time-reversal invariant systems whose classical analogues are K-systems 
show the same fluctuation properties as predicted by the GOE." 

^E-mail: nsw@ms.unimelb.edu.au 



Here GOE refers to the Gaussian orthogonal ensemble of random real symmetric matrices and K- 
systems are the most strongly mixing classical systems. The time reversal symmetry is responsible for 
the requirement that the random matrices have real entries. 

Another fascinating and enigmatic connection is that of random matrix theory and number theory, 
which concerns the statistics of the large zeros of the Riemann zeta function on the critical line. 
Incredible numerical evidence of Odlyzko ||^ based on the calculation of the critical zeros well into the 
asymptotic regime (zero number 10^° and 10^ of its neighbours) has convincingly shown their statistics 
to accurately follow the Gaussian unitary ensemble (GUE) of random Hermitian matrices predictions. 
This connection was first observed in analytic work of Montgomery Q on the pair correlations of the 
zeros. Denoting 

s„ = 1/2 + lE^ , d{E) = J26{E- En) , (1) 

n>l 

for the zeros and their density respectively, the asymptotic pair correlation was evaluated as 
1 T {d{E)d{E + e)) ^ sin^TTg 

= l + SGUE{i) (3) 

where e = e/{d{E)). This prediction, and other statistical measures of the asymptotic zeros such as 
spacing statistics, were tested in the numerical work of Odlyzko with almost perfect agreement being 
found at a graphical level. 

2 Matrix Ensembles, Statistical Mechanics of Levels 

We take as an example of a random matrix ensemble the GUE, which in application to chaotic quantum 
systems applies in the case of no time- reversal symmetry [H, T] 0. One constructs a. NxN Hermitian 
matrix X by taking real diagonal elements Xjj to be independently chosen with p.d.f. 

^e--« , (4) 

and upper triangular elements Xjk = Ujk + ivjk are independently chosen with p.d.f. 

TT TT 

The joint p.d.f. for elements of X is then 

N 

P(X) oc [] e-"-- n e-2|-^-l'= H e-l^^-l' =e-T^(^') . (6) 

j=l l<j<k<N l<j,k<N 

This distribution function is invariant under arbitrary unitary transformations 

P{U-^XU) ^ P{X) , (7) 

and hence the naming of the ensemble. The Gaussian form of the p.d.f. can be uniquely characterised 
by the invariance and the factorisation property 

P(.X) = l[fix,k) . (8) 

Alternatively it is the unique p.d.f. which maximises the entropy 

S[P] = - J ti{dX)P log P , (9) 

subject to the constraint (TrX^) = N^. Primary interest centres on the eigenvalue level statistics and 
so a description in terms of eigenvalues is required. This is achieved by the mapping 

, , , ^ I eigenvalues Ai < . . . < Aw , ^ 

i/27V(7V + 1) elements {x,, }-> J (10) 

I y2-™ (-™ — 1) variables Pi, ■ ■ ■ 
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where the pi are related to the eigenvectors. Using the transformation of Hermitian matrices to the 
diagonal representation the volume form is correspondingly transformed as 

N N 

{dX)= f\dxjj /\ dujk A dvjk ^ {UUU) Yl {Xj - Xkf /\dX, . (11) 

j = l l<j<k<N l<j<k<N 3 = 1 

This eigenvalue and eigenvector dependence thus factorises. The eigenvalue pdf can then be written as 

P(Ai,... ,AA,) = ^e-S.^' n (^.^ -^fc)'- (12) 

^ l<j<k<N 

This p.d.f. is proportional to the Boltzmann weight factor 

1 ^ 1 

from equilibrium statistical mechanics for a system of N charges on the line (or circle) at temperature 
1//3 located at positions Xj , and with the total potential energy 

U{xi, . . . ,xn) = ^ Inw(x-i) - ^ \n\xj-Xk\, (14) 

1<J<JV l<j<k<N 

(with — logw(a;) = x'^) composed of a 1-body confining potential (the first term) and a 2-body elec- 
trostatic interaction (the second term). Consequently a configuration integral can be defined 



oc 



Zn = dxi... dxNe-f^^ . (15) 



The gap probability, or probability for the exclusion of eigenvalues from an interval /, is a funda- 
mental statistic underlying the spacing distribution and is defined by 



I) = ^ fdxi... f dxNe-^"" , 

'^N J I J I 



(16) 



with / = (— oo, oo)\/, however this form is often not suitable to work with. Instead we seek a form 
involving the n-particle distribution function 



Pn{xi, . . . ,Xn) = ^ S{xi-Xj^)---S{Xn-Xjj'^ , 

dxn+1... dxjve-'^^f^^--^") . (17) 



A^! 1 



{N-n)\ Zn 

This can be done by introducing the generating functional 

N 



ZN[a] = — dxi... dxNT\[l + a{xi)]e-P''^^^'-'^-\ (18) 
Zn J-oo J-oo 

expanding out the product and using the definition ( p^ to obtain the formula 

ZAr[a] = l + ^— / dxi... dXnY\_0-{xi)Pn{xi,... ,Xn) ■ (19) 



1=1 

Choosing ~a{x) to be the characteristic function of / so that 



, , — 1 if a; G / , , 

«(^)Ho if x4l^ ^''^ 
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we see that ( Jig ) coincides with the definition (|16|) of Efj{0;I) and ( |19D then gives 

Ep{Q-I) = l + y^^-^ dxi... dXnPn{xi,... ,Xn) ■ (21) 
„=1 

The probabiUty -B/3(0; /) determines the p.d.f. that given there is a particle at oi, the nearest neighbour 
to the right is at 02, through the formula 

p^(0; (ai, a2)) ee f dx, . . . f dxN e-^^^^^-^^-^^- , (22) 

^ -^;3(0;(ai,a2)) • (23) 



p(ai) daida2 

The distribution p;3(0; (ai, 02)) is readily measured empirically, especially in the bulk region of the 
spectrum when translation invariance gives p^(0; (ai, 02)) = p/3(0; |ai — a2|). 

3 Fredholm Integral Operators 

For (3 — 2 the probability i?/3(0; /) can be written as the determinant of a Fredholm integral operator 
[^. To show this we make use of the fact that the general n-particle distribution has the determinant 
representation 

Pn{xi,... ,Xn) = det[ii:(xi,a;j)]"j^i , (24) 
where, with {pn{x)}n=f).i,... the orthonormal polynomials with respect to the weight function 'w{x), 

K{x, y) = [w{x)w{y)Y/^ ^ Pn{x)pn{y) , (25) 
With a„ denoting the coefficient of x" in Pn{x) the Christoffel-Darboux summation gives 

j^. X OJV-lr f X f if2 PN{x)pN-l{y) - PN-l{x)pN{y) 

K{x,y) = [w{x)w{y)\ ' , 



aN x~y 

Hx)i^{y) - (f>{y)i'{x) 
x-y 



(26) 



The determinant formula (E4) substituted in (21) gives a well known [p| expansion formula for the 



determinant of an integral operator in terms of its kernel. Explicitly 

E{0; I) = det(I - K) (27) 

where IK is the integral operator on / with kernel K(x, y). 

Tracy and Widom 0| have shown that by introducing a number of auxiliary quantities related to 
K, for some special choices of w{x) and / the Fredholm determinant ( ^ ) can be specified in terms of 
the solution of a nonlinear equation. The auxiliary quantities required are the operators 

= R(x,y) , = p(x,y) , (28) 

where the notation A = A{x,y) denotes that the integral operator A has kernel A{x,y); the pair of 
functions 



Qix) = / dy p{x, y)4>{y) = (I - K)- V , 

^\ (29) 
P{x) = / dy p{x, y)ij{y) = (I - K)" V , 



qj = lim Q{x) , pj = lim P{x) ; (30) 
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and the scalar products 

u = {^\Q) - / dy Qivmv) 



I 

V = (^IQ) = j^dy QivMy) = = J^dy P{y)c^{y) , (31) 



w = (V^IP) = / dy P{y)iAy) 



I 

The essence of the Tracy & Widom machinery in this setting is to set up a system of partial 
differential equations for the above quantities with respect to the endpoints of the interval, defined by 
/ = (oi, 02). One set of equations are universal while a second set of differential equations describes the 
coincident functions and their derivatives, and includes specific features of the particular orthogonal 
polynomial system at hand. In particular one requires the differential-recurrence relations 

m{x)(f)'{x) = A{x)(j){x) + B{x)^p{x) , 

m{x)ij'{x) = ^C{x)(t){x) - A{x)il){x) , ^ ' 

where m, A, _B, C are polynomials in x parameterised by 

■m{x) = fia + t^ix + ii2X^ , 
A(x) = ao + aix , 
B{x) ^ /3o + Pix , 
= 70 + lix , 

which is a common structure for the classical orthogonal polynomials. 



4 Spacing p.d.f. and Painleve transcendents 

The classical weights are 

{e~^^ , Gaussian 
x°-e~^, X > 0, Laguerre (34) 

(l-a;)'^(l + a;)^ -Kx < 1, Jacobi, 

which together with their corresponding orthonormal polynomials give rise to functions (j){x), ip{x) sat- 
isfying ( p2| ) and (^3|). For general parameters in all these cases, it has only been possible to characterise 
£'2(0;/) in terms of a nonlinear equation when the interval / includes an endpoint of the support of 
w{x) and so depends only on one variable. It turns out that then the auxiliary quantities q' /q, p' /p and 
R = R{s, s), as well as the probability £'2(6;/) can be expressed in terms of Painleve transcendents. 

In Tables 1-3 we list the evaluations for the three classical weights respectively. Also listed are 
the specific coefficients (p3|), and an equation or equations (integrals of motion) relating the quantities. 

To illustrate the construction of these tables, consider the Jacobi weight with / = (—1, s). Direct 



application of the Tracy and Widom theory gives the coupled differential equations 

[ln£;2]' = , (35) 

(1-5^)9' = [ao + ais + v]q+ [Po + u{2ai-l)]p , (36) 

(l-s^)p' = -[70 - w{2ai + l)]q - [uq + ais + v]p , (37) 

u'^q\ (38) 

v' =qp , (39) 

w' =p^ , (40) 

{l-s^)R = [70 - w{2ai + l)]q^ + [(3o + u{2ai-~l)]p^ + [ao + ais + v]2qp (41) 

[{l~s^)R\' ^2aiqp , (42) 



From this system of differential equations two integrals can be constructed , which are given in Table 
^, and using these a sequence of eliminations leads to a second degree second order differential equation 
for cr(s) = {l — s^)R{s). This was first found in an indirect manner in In Table |^ we identify the 
solution of this equation with a number of P-VI transcendents, and give that specific mapping in the 
first case listed. For this purpose we use theory from Q]. 
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5 Conclusions 



In this work we have considered the evaluation of the gap probabihty i?2(0; /) for the classical weights 
with / including an endpoint of the support. This quantity, and the auxiliary quantities ll 
and R = R(s,s), are given in terms of Painleve transcendents thereby extending findings of earlier 
works 1^, ^. Our approach has been to construct all the integrals of the motion which allows us 
to reduce the system of differential equations arising in the Tracy and Widom formalism to a single 
second order ODE by the most direct path possible. Furthermore we have displayed the appearance of 
the Painleve transcendents in the random matrix context in a transparent manner, which we consider 
to be essential for an understanding of the relevance of integrability in this situation. 
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Gaussian weights and Hermite orthogonal polynomials 
w{s) = e"^', m{s) = 1, Ais) = -1, B{s) = C{s) = V2N 
I = (s, oo) 

V2N[u — w] + 2uw = qp 



E2{s)=exp{-J^dtR{t)) 



R{s) 



■)N-1 



q' N 

— ~ — s H 

q s— >oo s 

p' N-1 

— - -s+ 

p s^oo s 



u;{s) PIV a = 2N -1, /3 = 



Q 

1 = 
P 

R = 



2Nuj 



= —S — LJ — 



-1/20; + -c' 



-1/2(3' -2N)lO - i/2SL0' - i/sw' + — (w') 



l\2 



w(s) PIV a = -N±l, (3=-2N^ 



R 



TV^ 1 



Table 1: Painleve transcendents in the Gaussian Unitary Ensemble 
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Laguerre weights and orthogonal polynomials 

w{s) = s^e-' 
m{s) = s 

A{s) = -1/20 -N + 1/2S 
B{s) = C{s) = ^JN{N+a) 



/=(0,s) 



\/ N{N +a)[w — u]+ WW = sqp — sR 



E2{s) = exp {- dt R{t)) 



a{s) = sR{s) 



T{N + a+l) 



(Af-l)!r(a+l)r(a + 2) 
q' a 2N-l+a 
q s-yQ 2s 



a+1 - 

s e 



2N-2 

1 TTS 

a+2 



P' 



2{a+l) 
27V+l+a 



p s^o 2s 2(a+l) 



uj{s) PV 



a 




13 = 


7 = 


6 = 


1/2(1 - 


-af 





-2N-a 


-V2 




V2 




2N+a 


-V2 


1/2(1 -a- 


Nf 


-1/2^' 


—a 


-V2 


y2(l- 


Nf 




a 


-V2 



q'_{a-l)u)-a 2Nlj u'-l 

'q ~ 2^ {lj+lj')s + (a-l)w(a;-l) 2{uj~l) 

p'_{a-l)uj-a 2{N+a)uj uj' + l 

p" ^ 2^ ^ ' 



sR= - — 



su 



{uj-Lo')s - {a-l)uj{uj-l) 2(w-l) 

n 2 



— U 



+ y^a^uj + y.2{2N+a)s—^ + 1/45 



(a;- 1)2 



Table 2: Painleve transcendents in the Laguerre Unitary Ensemble 
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Jacobi weights and orthogonal polynomials 

wis) = , m{s) = l-s\ A{s) = ^^^^"^^^^ - ^^+"+^ = ao + ais 

, . ^ 2^N{N+a){N+b){N+a+b) fW+a+b+l 

2N+a + b \l 2N + a+b-l~ ^° 



_ 2y/N{N+a){N+b){N + a+F) j 2N + a+b-i ^ 

2N+a+b \l 2N+a+b+l^^° 



I=i~l,s) 



cr(s) = (1 - s'^)R = -{2N+a+b)v 
[/3o + u(2ai-l)][7o - w;(2ai + l)] = /3o7o - (l-s^)o-' - sa + —a + -^cr^ 

Of ]^ 40!^ 



r{N + a+b+l)r{N+b+l) fc^i / 2N^+2N{a + b)+ab-b 

''^^> sS^-i 2b{N-iy.riN+a)T{b+i)rib+2y'' ' { 2(6+2) ' 

7-2 



b _ 2N^+2N{a+b+l)+a{b+l) 
7 2(s+l) 4(6+1) 

y 6 _ 2Af^ + 2Ar(a+6-l)+a(6-l)-26 

p" s^-i 2(s+l) 4(6+1) 



a;(a;) PVI a = 

V2 

i/2[l - (2iV+a+6)]2 

i/2[l-Af-a-6]2 
i/2[l-iV]2 
y2[l-A^-a]2 
i/^[l-N-bf 



7 = 



-V2«' 1/26 



i/2(2iV+a+6)2 
-i/2(27V + a+6)2 

W + bf 

-y^iN+a+bf y^N^ 



-l/2(7V + 6)2 



6 = 

y2[l - {2N+a+bf 



y2 

V2[l - a^] 



i/2[l - N'] 
i/2[l - (7V+a+6)2] 

y2[i - {N+bf] 

V2[l - (iV+a)^] 



s = 2x-l 

(l-s^)^ = + ^ + (27V+l+a+6)^^^^ + x{l-x)- 



w — a; 



a; - 1 + w - (27V - 1 + a + 6) ^^^-^ + a;( 1 - x) 



u — x 

' U) — X 



a = 



2lij{liJ—1){u! — x) 



x{x—l) 



1/20'- + 1/26' 



X— 1 
LU-1 



y2(27V+a+6) 



2X(1-X) 



Table 3: Painleve transcendents in the Jacobi Unitary Ensemble 
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